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Abstract 

The phase structure of hadrons and quark-gluon plasma is investigated by two 
types of equation of hadron state, namely ideal hadron gas model and the com- 
pressible bag model. It is pointed out that, while the ideal gas model produces 
unrealistic extra hadron phase, the compressible bag model gives an expected and 
reasonable phase diagram even if rich hadron spectrum is taken into account. 
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1 Introduction 



The quark-gluon plasma (QGP) has been expected to appear at high temperature and/or 
high density, due to asymptotic freedom of QCD. Indeed, not a few of the experimental 
results upto CERN SPS incident energies suggest the existence of the QGP phase and 
the recent experiment at BNL RHIC may give us cleaner and richer signals ||. 

Theoretically, in the lattice QCD calculations @, £|, |5j, || [7j, no phase transition is observed 
for the two light or two light and one medium light flavor cases at zero baryon density. In 
both cases only a continuous cross over is observed for realistic quark mass values. 

On the other hand, in more phenomenological approach, state equations of QGP and 
hadrons are assumed, and the transition point is determined by Gibbs condition. In this 
approach, the cross over cannot be reproduced, and that may give rise to a reason to doubt 
its validity. However, the authors believe that the gross structure of the phase diagram 
obtained in this approach should still remain valid far away from the cross over region. It 
should also be noted that at a finite baryon chemical potential, not much is known from 
lattice QCD because of a well known technical difficulty which makes the Monte Carlo 
technique inapplicable. 

The purpose of the present paper is to investigate what type of the state equation of 
hadrons is preferable in the phenomenological approach, where the discussion is devoted 
exclusively to gross phase structure of hadrons and QGP, neglecting the multi-quark states. 
It means that the detail structure around phase transition or cross over region is neglected. 

In the phenomenological approach, there are following problems concerning the gross 
structure of phase diagrams. When a simple ideal gas of nucleons is used for the equation 
of hadron state and MIT bag model is used for the equation of QGP state, hadron phase 
appears at high density and at zero temperature (see, e.g. ||). This is because 

pE/p;- 27 as Hb^og, (T = 0) (1.1) 

where p° h and p° denote pressures of ideal nucleon phase and (ideal) quark phase, respectively, 
and hb is baryochemical potential. To circumvent this difficulty, volume exclusion effects 
of hadrons of van der Waals type were taken into account (see, e.g. |8|, see also |T0[] ). 
However, the treatments in || |9|, |K| are thermodynamically inconsistent fill , for example, 



n 7^ (dp/dfi)T, s 7^ (dp/dT)^. Although the formulation in |TT| is thermodynamically 



consistent, it has still a difficulty that pressure diverges as n — > 1/vq, where n is number 
density and vq is fixed hadron volume, and acausality emerges (there is no relativistic rigid 
body) as the authors themselves pointed out. Thus we were led to 'soft core' model, that is 
the compressible bag model [O, O. 
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In addition to the difficulty that hadron phase appears at high density and at zero 
temperature, an equation of simple ideal gas of hadrons again suffers from the same difficulty 
at another region, that is, at high temperature and zero baryon number density jn] . The 



argument in [13] is simple. At /is = and in high temperature region, the pressure p q of 



QGP phase and the pressure ph of hadron phase are given by 

p q = ^(l/90vr 2 )T 4 -B, p h = g h (l/90n 2 )T\ (1.2) 

where g q and are the degeneracy factors in each phase and B is a bag constant. When 
we consider a mixed gas of lower lying nonstrange hadrons, that is, n, p, N and A, the 
degeneracy factor g^ becomes 3+9+(7/4)(4+16) = 47 and it exceeds g q = 2-8+(7/4)4-3 = 37. 
Thus hadron phase appears again. In the compressible bag model, however, the effective 
degeneracy factor of hadron phase is greatly reduced so that QGP phase is realized at high 



temperature even if infinitely many nonstrange mesons are taken into account [15| 



The authors of [|T|, [13], [H| suggested that the compressible bag model gives the phase 



diagram free from above difficulties, but it has not shown explicitly: In [T3| the equations 



of state only at T = are used for discussing neutron stars, and in []L5] used are those at 



Hb — for multiple production in heavy ion collisions. In [y| the phase diagram with finite 
T and /xg is calculated, but rich hadron spectrum is neglected. 

The unique purpose of the present paper is, then, to show that the compressible bag 
model indeed gives the naively expected phase diagram in all regions. As far as the authors 
know, it is the only model that satisfies the followings: 

(a) It is built in thermodynamically consistent formalism. 

(b) It is valid in very high temperature regions even if rich hadron spectrum is taken into 

account. 

(c) It is explicitly calculable in whole regions, and it gives well-behaved phase diagrams, i.e. 

there are no extra hadron phase. 

Since the hadron level structure affects the phase diagram, we consider the following two 
models: 

model I: Hadron phase consists of NNtt system while QGP phase consists of nonstrange 
quarks, their anti-quarks and gluons. 

model II: Hadron phase consists of a system of 102 species of hadrons of which masses are 
less than 10.6 GeV, except for hadrons with top quark |T(J. QGP phase consists of u-, 
d-, s-, c-, b-quarks and their anti-quarks and gluons. 
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Model I is a model with a few hadrons, while model II is an example which includes many 
hadrons. 



While in |T5[ the continuous level-density function of [T7| is used, in this paper the 



phase diagrams are calculated with real spectrum of hadrons in order to avoid the model 
dependence. The effect of the cutoff in the hadron spectrum in model II is discussed in the 
final section. 

In section 2, for comparison the phase structure in free point-like model is shown with 
the above two types of models, where one can explicitly see the difficulties explained in this 
section. In section 3, the equation of state in the compressible bag model [|13|, [L5| is briefly 
reviewed, and with the resulting phase diagram it is shown that the difficulty is indeed 
removed. Final section is devoted to concluding remarks. 



2 A difficulty in the phase structure for models of free 
point-like hadrons 

In this section, we assume that all hadrons, quarks, and gluons are free point-like objects 
and discuss their phase structure. For this purpose, we present expressions of pressures of 
free point-like particles in order to establish notations in the first subsection. On the basis 
of these expressions, the phase diagram is discussed and a difficulty is pointed out in the 
second subsection. 



2.1 Pressures for hadron gas and QGP 

Total pressure p of a mixed gas of iV-species of free point-like particles is given by 



N 



(2.1; 



n=l 



where T (p, B ) is the temperature (baryochemical potential) of the gas and p n is the 
pressure (baryochemical potential) of the n-th species of particles. The pressure p n is given 
by 

d 3 k 



Pr, 



Vng n T 



(27T) 



log 



1 + r] n exp 



f 



k 2 + m 2 , 



(2.2) 

(2.3) 
(2.4) 
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where rj n is a statistical factor that takes the value —1 for bosons and +1 for fermions, 
respectively. The quantities g n , a n , and m n denote the degeneracy factor, the baryon number, 
and the mass of the n-th species of particles, respectively. 

In low temperature region, p n is expanded as follows; 

/in + k n f 



Pn 



9n 



24tt 2 



9n , rr.2 , 7T[2 9 



^m 2 n fi n k nf + \fi n log 



k 



nf 



+—fi n k nf T 



720 



3/i n 

k n f 



knf, 



m r , 



(2.5) 
(2.6) 



when the n-th species of particles are fermions. In the above k n f is a fermi momentum. 
Eq. (|2.5|) is applicable for fi n > m n . Otherwise the pressure p n vanishes identically. 



Now let us present expressions for pressures of hadrons and QGP. As for total pressure 
Ph of a mixed gas of iV-species of free point-like hadrons, it is given by 

N 

p h {T,fi B ) = Y,PhA T >Vi)> (2-7) 

i=l 

where follows (|2.2|) ~ (|2l]|) . In the followings we consider two cases for hadronic systems; 
that is, the model I and the model II which are defined in the introduction. 



On the other hand, total pressure of QGP p q is given by 



N' 



p q {T,ii B ) = J2p qj (T^j) - B, 



(2.1 



in the bag model. Here p q . follows (2~2)~(2!d) and B is a bag constant. The index j denotes 
the species of particles, and we consider two cases, the model I and the model II. 



2.2 Phase structure 

Given the expressions of pressures of hadron phase and QGP phase, we can consider their 
phase structure. Critical points are determined by the following Gibbs condition; 

Ph(T,fi B ) =p q (T,fi B ). (2.9) 

In Fig. 1, ph and p q are plotted as functions of fi%, where T is fixed to zero. In Fig. 2, they 
are plotted as functions of T 4 , where \i B is fixed to zero. In both figures, calculations are 
performed for the model I. As shown in Fig. 1, there is no critical point for B 1 ^ > 0.1536 GeV. 
For = 0.1536 GeV, there is one critical point. For 0.1534 GeV < B 1 ^ < 0.1536 GeV, 
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Fig. 1: p h and p q as functions of /ifg for the 
model I (T = 0). The unit of 5 1 / 4 is GeV. 



Fig. 2: and p q as functions of T 4 for the 
model I (fi B = 0). The unit of B 1 ' 4 is GeV. 



there are two critical points. For I? 1 / 4 < 0.1534 GeV, there is one critical point since ph 
vanishes for fi B < m N and unphysical is an apparent critical point in the region \i B < mjv, 
where is nucleon mass. Therefore, at T = 0, abnormal hadron phase always appears at 
high densities. 

The reason for this difficulty is as follows. In higher density region of fi B 3> m Ni masses 
are negligible so that pressures ph and p q are approximated as 



Ph ~ 9hPB> ( 2 - 10 ) 
V q ~ 9 q f4, (2.H) 

where (g g ) is a statistical degree of freedom of hadrons (QGP). Since \i q = (1/3) fi B , ( |2.11[) 
reads 

p q « (2-12) 

Then becomes large faster than p q &s fi B increases, since an effective statistical degree of 
freedom of QGP is greatly reduced. In other words, hadrons get a share of energy three times 
larger than QGP, when fi B increases, so that ph becomes larger than p q . As a consequence 
hadron phase appears in high density region. 

Since the difficulty mentioned above stems from the qualitative nature of the model of 
hadrons that their statistical degree of freedom is too large, we should consider to modify 
the model of hadrons. So far we consider the model I at present, the situation will become 
worse in the model II, since the statistical degree of freedom becomes further large. 

In Fig. 3, phS for the model I and the model II are plotted as functions of /j, b , where T 
is fixed to zero and in Fig. 4, p h 's for those are plotted as functions of T, where \i B is fixed 



5 



0.8 



0.6 



> 
u 
O 

< 



0.4 



0.2 



20 



model I 
model II 



0.5 1 1.5 

H B [GeV] 



15 



> 
u 

O 



ft, 



model I ■ ■ 




model II 


/ 

/ 
/ 








/ y 

/ s' 






i 



0.05 



0.1 

T [GeV] 



0.15 



0.2 



Fig. 3: ph as functions of fi for the model I Fig. 4: ph as functions of T for the model I 
and the model II (T = 0) and the model II = 0) 



to zero. As seen in Fig. 3 (Fig. 4), p^ for the model II becomes considerably larger than p^ 
for the model I in high density (temperature) region, as expected. 

In Fig. 5, ph and p q are plotted as functions of fi% with T = 0. In Fig. 6, they are plotted 
as functions of T 4 with fig = 0. In both figures, calculations are performed for the model 
II. Qualitative features in Fig. 1 and Fig. 5 are the same. However, qualitative features in 
Fig. 2 and Fig. 6 are different. As shown in Fig. 6, hadron phase appears in high temperature 
region. The resulting phase diagram for the model I is shown in Fig. 7 and that for the 
model II is shown in Fig. 8. The situation for the model II has become worse than for the 
model I. In order to modify this situation, we should replace the model of hadrons with 
alternative ones. As one of them, we formulate the compressible bag model in the next 
section. 



3 The compressible bag model 



First, we briefly summarize the results of |13|, Let us suppose a mixed gas of n-species 
of hadrons enclosed in volume V at finite temperature T. In the compressible bag model, 
the free energy function F of the gas is given by 

F = F f (Ni, V, T, Mifa)), (3.1) 

i=l 

n 

V' = V-b^NjVj, (3.2) 
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where Ni is a number of z-th species of hadrons and V{ is its volume and Mj its mass. The 
constant b is a volume exclusion efficiency parameter. The function Ff is a free energy 
function of free point-like hadron gas. As for mass function M{(v i), we assume that of MIT 
bag model: 

Mi(vi) = A iV ~ 1/3 + Bv u (3.3) 

where B is bag constant. 

Under the approximation that the average of inverse Lorentz factor of z-th species of 
hadrons equals to unity ~ 1) in the rest frame of the system, basic requirement of the 

compressible bag model that dF/dvi = and the requirement that the chemical potential of 
hadron should be, /ij = dF/dNi = a^B, where a, is the baryon number of the i-th species of 
hadrons, determine the pressure p of the system as a function of T and /ie by the following 
equation: 

d 3 k 



p 



n r d 3 k 

VidiT J j^\og{l +r ]i eM-{Ei-^{p,rn i ))/T]}, (3.4) 



E l = Jtf + M^m,) 2 , 



M t (p, mi ) = ™*(l + ff) + 3/ , (3-5) 

bjn>ijp ( bp \ 

Viip, m i) = IH- bviV = aiHB ~ -j^- I 1 + -g 1 , (3.6) 

where gi is a degeneracy factor of the i-th species of hadrons and m; = 4(Aj/3) 3//4 5 1//4 is its 
mass in the vacuum. 

Here two comments are in order. First p determined as above depends on b and B only 
in the combination oib/B since Mj and \J { depend on b and B in that combination as seen 
in ( [3.5D and (|3~6| ). Second the pressure in the compressible bag model does not become so 
large in high density and/or high temperature region in contrast to the one in the model 
of free point-like hadrons. This is because masses of hadrons become large in high density 
and/or high temperature region in the compressible bag model so that kinetic energies of 
hadrons increase rather slowly and the pressure of hadrons do not become so large. Thus 
the compressible bag model have a chance to evade the difficulty discussed in the previous 
section. In the following, we assure that the difficulty is indeed removed, by numerical 
calculations. 

In order to do numerical calculations, we have to fix bag constant B and volume exclusion 
parameter b. The parameter b is determined by the relation 

b»N(p = 0) = ^j- = f (0.82 fm) 3 , (3.7) 
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Fig. 9: (a) A solid line shows ph as a func- 
tion of hb for the compressible bag model. 

(b) A dashed line shows p q for QGP phase. 
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and T = 0. 
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when the bag constant is given. The value 0.82 fm is the proton charge radius. It should be 
noted that equations of state for hadron phase are determined uniquely at this stage since 
the equation of states for hadron phase only depend on b/B. We tentatively put 

B 1 ' 4 = 0.25 GeV (3.8) 

and then get b = 5.0. The resulting ph and p q are shown as functions of \ib (T) in Fig. 9 
(Fig. 10). We can recognize that the compressible bag model does not suffer from the 
difficulty that exist for free point-like hadrons. This point is more clearly seen by the phase 
diagrams shown in Fig. 11 with B 1 ^ = 0.25 GeV. Critical temperature T c and critical 
baryochemical potential /i c are given by 

T c = 0.18 GeV, /i c = 2.2GeV in model I, (3.9) 
T c = 0.17GeV, /i c = 1.8GeV in model II. (3.10) 

It is noted that the state equation of hadrons are determined by the ratio b/B and that 
of QGP by B, the critical values, then, depends on B and b. If the constraint ( |3.7| ) is taken, 
the values depend only on B. As a reference, some cases are shown below; 

T c = 0.14GeV, /i c = 1.7GeV in model I, (3.11) 
T c = 0.14 GeV, fi e = 1.4 GeV in model II, (3.12) 



for B 1 ^ = 0.20 GeV and 



T c = 0.22GeV, /i c = 2.7GeV in model I, (3.13) 
T c = 0.21 GeV, /i c = 2.1GeV in model II, (3.14) 



for 5 1 / 4 = 0.30 GeV. 



From Figs. 9~11 one can see that the compressible bag model gives an expected and 
reasonable phase diagram in whole regions. This result is stable in the sense that one can 
choose the bag constant from rather wide range. 



4 Concluding remarks 

In this paper, it has been shown by explicit calculations that the compressible bag model 
gives well-behaved phase diagram in whole regions, even if many hadron states are taken into 
account. In the present calculation it is true that, even in the model II, the infinite series 
of hadrons are cut off at finite mass, but the compressible bag model may give an expected 
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phase diagram if infinite series of hadrons are included. This is because, in the model, the 
masses of hadrons increase and the effect of higher mass state is much suppressed in high- 
temperature or high-density regions. Indeed in a continuous level-density function are 



used for hadrons and it is shown that the abnormal hadron phase does not appear in the 
compressible bag model, although limited at \i B = 0. 

Free point-like models produce unnatural hadron state in high-temperature or high- 
density regions as shown in Sections 1 and 2. A way to avoid this difficulty is the compressible 
bag model, and it is worth examining in various phenomenological analyses. 
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